For a graph G=(V,E) with v(G) vertices the partition function of the random cluster model is defined by ZG(q,w)=∑A⊆E(G)qk(A)w|A|,\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\begin{aligned} Z_G(q,w)=\\sum _{A\\subseteq E(G)}q^{k(A)}w^{|A|}, \\end{aligned}$$\\end{document}where k(A) denotes the number of connected components of the graph (V, A). Furthermore, let g(G) denote the girth of the graph G, that is, the length of the shortest cycle. In this paper we show that if (G_n)_n is a sequence of d-regular graphs such that the girth g(G_n)rightarrow infty , then the limit limn→∞1v(Gn)lnZGn(q,w)=lnΦd,q,w\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\begin{aligned} \\lim _{n\\rightarrow \\infty } \\frac{1}{v(G_n)}\\ln Z_{G_n}(q,w)=\\ln \\Phi _{d,q,w} \\end{aligned}$$\\end{document}exists if qge 2 and wge 0. The quantity Phi _{d,q,w} can be computed as follows. Let Φd,q,w(t):=1+wqcos(t)+(q-1)wqsin(t)d+(q-1)1+wqcos(t)-wq(q-1)sin(t)d,\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\begin{aligned} \\Phi _{d,q,w}(t):= & {} \\left( \\sqrt{1+\\frac{w}{q}}\\cos (t)+\\sqrt{\\frac{(q-1)w}{q}}\\sin (t)\\right) ^{d}\\\\{} & {} +\\,(q-1)\\left( \\sqrt{1+\\frac{w}{q}}\\cos (t) -\\sqrt{\\frac{w}{q(q-1)}}\\sin (t)\\right) ^{d}, \\end{aligned}$$\\end{document}then Φd,q,w:=maxt∈[-π,π]Φd,q,w(t),\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\begin{aligned} \\Phi _{d,q,w}:=\\max _{t\\in [-\\pi ,\\pi ]}\\Phi _{d,q,w}(t), \\end{aligned}$$\\end{document}The same conclusion holds true for a sequence of random d-regular graphs with probability one. Our result extends the work of Dembo, Montanari, Sly and Sun for the Potts model (integer q), and we prove a conjecture of Helmuth, Jenssen and Perkins about the phase transition of the random cluster model with fixed q.